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Abstract 

The form of an effective electron-electron interaction in a quantum wire 
with a large static dielectric constant is determined and the resulting prop- 
erties of the electron liquid in such a one-dimensional system are described. 
The exchange and correlation energies are evaluated and a possibility of a 
paramagnetic-ferromagnetic phase transition in the ground state of such a 

system is discussed. Low-energy excitations are briefly described. 
Pacs: 71.10.-w, 71.10.Pm, 73.23.-b 
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In three and two dimensional systems any interaction between electrons (in a forward 
scattering channel) appears to have a marginal property.0 This means that the low-energy 
dynamics of such systems is described by the Fermi liquid theory where quasi-particles 
are in a one-to-one correspondence with bare electrons^ One-dimensional systems are very 
exceptional in this respect. The interaction between fermions is relevant and leads to a 
new collective type of a particle motion at low-energies. Such a system is described by 
the Luttinger liquid theory, according to which a spin-charge separation and an anomalous 
scaling take place.! 

A recent progress in nanofabrication enables to construct various one-dimensional struc- 
tures, in which the Luttinger liquid theory can be experimentally tested.0 In such structures 
the effective interaction between the electrons is not a short-range one as in the standard 
Tomonaga-Luttinger (TL) model.! Rather it must be of a long-range type because it origi- 
nates from the Coulomb electrostatic forces between the charged particles. Accordingly, in 
such a case certain properties of the Luttinger liquid will be modified. 

In this Brief Report we analyze theoretically the form of an effective interaction between 
the electrons in a one-mode quantum wire, and the corresponding new properties of the 
system. In particular, we are interested in how ground-state and excitation energies will be 
affected by a difference between the dielectric constants of the wire and its environment. Such 
a difference is large in the case of free standing wires, particularly of IV- VI semiconductors 
such as PbTe, characterized by a large static dielectric constant e ~ 1000.1 This very large 
dielectric constant is due to the proximity of this system to the ferroelectric phase transition. 
In those cases the electrostatic potential between the electrons is strongly modified because 
of the presence of image charges that assure the correct boundary conditions. 

In order to determine the form of the interaction in such a system we consider a model, 
in which the electrons can propagate along an infinitely long cylinder of radius a, made of 
a material with a macroscopic dielectric constant e\. This cylinder is embedded into a bulk 
system with a macroscopic dielectric constant €2- Further, we assume that the single-electron 
wave function vanishes at the boundary of this cylinder. In other words, an electron moves 
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in a cylindrical potential well with infinite barriers. Eigenfunctions ^ n mk in such a geometry 
are readily to find, and the result in the cylindrical coordinates (p, 4>, z) is, 
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where J m {x) is the Bessel function,! X mn are nodes of J m (x), (n,m,q) is a set of quantum 
numbers. 

In the following we are interested in properties of the quantum wire when only the lowest 
quantum level (n = m = 0) is occupied. In this case the radial dependence of the wave 
function (1) can be approximated very well by the following parabola! 

^(1 - (£) 2 ) for p<a 
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otherwise. 
An electrostatic potential V(r — r ) between two electrons is determined by solving a 
Poisson equation V 2 V = ^|p<5(r — ro), where e is an charge of an electron, with appropriate 
boundary conditions at p — a, i.e., V\ = V2 and eid p V\ = e2<9 p V2.! As a result we find that 

e 2 00 



Hp) 



Y, / d\e im ^- M cos[\(z - z )} 



(X P< )K m (\p > ) + 



'I — q\ w 



K' m (\a)K m {Xa) 



-Im(Xp) 



fK' m {Xa)I m {Xa)-K m {Xa)I' m {Xa^' n "'' r> ' (3) 
where p< = min(p, p ), p> = max(p, p ), I m (x), K m (x) are the modified Bessel functions,! 
and primes denote their derivatives. We have not considered here effects coming from 
electrodes, which are usually attached to the system, assuming that they are very far from 
each other and modify the system properties only very close to the edges. However, for short 
wires and quantum point contacts the electrodes may be important as well. 

An effective Hamiltonian describing electrons in a ID quantum wire with the lowest level 
occupied only (n — m — 0) has a two-body matrix element (determined in a one-particle 
basis of states (2)) in the following form: 
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where 
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Mlaql) = ^(M)«-^(hl)«)' (5) 

The first part, which was found previously in Ref. 7, corresponds to the long range Coulomb 
interaction between the electrons moving in a quasi- ID constriction. The second part, which 
disapears when e\ = e 2 , describes the interaction between an electron and the image charges, 
which assures the proper boundary conditions. 

In Fig. 1, we plotted the matrix element V(aq) as a function of aq for different ratios 
£1/^2 = 1, 10, 100. We see that V(aq) is a decreasing function of aq and for large values of aq 
becomes relatively small. Additionally, we find that the numerical values of V(aq) diminish 
when the dielectric constant e\ increases and this reduction is very different when the image 
charges are taken into account as is shown in the inset to Fig. 1. For small aq, however, the 
effective interaction is seen to diverge. We have been able to examine this limit analytically 
and found that at \aq\ <0 the matrix element behaves as 
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so it diverges logarithmically (7 is the Euler constant). This behavior is characteristic for a 
ID Fourier transform of the Coulomb interaction e 2 /r (with a as a short-distance cut off). Bi 
However, it is surprising that V(aq) for \aq\ — > does not depend on the dielectric constant 
ei of the wire but only on the dielectric constant e 2 of the environment. This means that as 
the wire is infinitely thin (or particles are very far from each other) a single electron interacts 
mainly with the image charges and not directly with the other electrons in the wire. We 
have checked explicitly that the disapearence of the dielectric constant ei in the first leading 
term takes place for both wave functions (|1|) and (0), which are separable in the cylindrical 
coordinates.© 

Having determined the analytical form of the electron-electron interaction matrix element 
we can calculate ground state properties of such a quantum wire. We assume that the 



electrons propagate in a jellum environment with positive ion charges. This leads to the 
cancelation of a direct self-energy contribution, which otherwise brings about infinities in 
the perturbation expansion.^ Hence, the Hartree-Fock correction to the ground state energy 
has the simple form 

and its numerical values are depicted in Figs. 2a and 2b, for the different radii (a) of the wire 
and the different RPA parameters r s = (2aQn)~ 1 , where n is the density of the ID electron 
gas.0 All results below are presented in the atomic units where = is the effective 
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Bohr radius, and Ry* = 2 rrf _ is the effective energy unit; so-called the effective Rydberg. 
Note that since in the leading term of the expansion only the dielectric constant e<i 
appears, we have defined the effective Rydberg and the effective Bohr radius with e^. We 
see from Figs. 2 that, as expected, the absolute value of the exchange energy is smaller for 
either a wider wire (a = 10), a less dense electron gas in the wire (r s = 4), and a greater 
dielectric constant e\. However, the decay of Efj F with e\ is significantly weaker in our model 
(f|) than in the model of Ref. 7, in which no image charges were considered, A = 0, as is 
shown in the inset to Fig. 2a. This comaprision convincingly demonstrates the important 
influence of the image charges upon the ground state energy of the wire. 

We have also evaluated a correlation energy E c , i.e., the correction to the Hartree- 
Fock energy originating from a linear screening of the electron-electron interaction due to 
electron-hole excitations. For that we calculate the effective (RPA) interaction, which in 
the static limit has the form0 V^J A {aq) = V(aq)/e(aq), with the dielectric function e(aq) 
approximated as e(aq) ~ 1 + 2V{aq)/ r KVp, where vf = is the Fermi velocity. The 
correlation energies E c are plotted in Figs. 3a and 3b for the various wire's radii a and the 
parameters r s . 

In order to estimate the transition line between the paramagnetic and the ferromagnetic 
ground states we compare the ground state energies for these two phasesiil Our results are 
shown in Fig. 4 for the two wire's radii a = 1 and 10. The critical value of the RPA parameter 

5 



Tg, above which the system would be completely polarized, increases with increasing e\. 
This increase is, however, sublinear which indicates that the effects of the Coulomb carrier- 
carrier interaction may remain important even in nanostructures with a large bulk dielectric 
constant. At the same time, if there was e\ <C €2 then the critical value of r c s would go 
rapidly to zero (c. f. the inset to Fig. 4). This result suggests that if a quantum wire with a 
small dielectric constant was deposited on a ferroelectric substrate, then the ferromagnetic 
phase instability would be even more likely. 

It is tempting at this point to make a comment on the so-called 0.7 step in the quantized 
conductance, which is observed in the quantum point contacts of GaAs/AlGaAs,E! and 
also in PbTe.i Our results make possible to evaluate electron concentrations at which the 
zero-temperature ferromagnetic instability might appear. However, we cannot exclude other 
types of instabilities, e.g., the charge- or the spin-densities- waves, which are not discussed 
in the present paper. 

Finally, we briefly describe the low energy excitations in our model. As is known, ar- 
bitrary weak interaction destroys the Fermi liquid description in one dimensional systems. 
Instead, the Luttinger liquid theory emerges as the proper low-energy principle in this case.H 

The one-dimensional model of interacting electrons can be solved exactly in the low- 
energy limit by means of a bosonizationi. In our case we must linearized the dispersion 
relation — ep = vp(±k — hp), where vp = hkp/m* is the Fermi velocity at the two Fermi 
points ±kp. Next, we introduce the so-called left and right moving operators corresponding 
to ±vp, respectively, and then define fluctuation density operators p qaa = J2k c k+qaa c kaa for 
each branch a = R,L separately. In terms of these operators the many body Hamiltonian is 
bilinear and can be diagonalized exactly because in the low-energy limit p qaa obeys boson- 
like commutation relation.! 

As a result we find that spectra in the charge and the spin channels are different . Namely, 



the eigenvalue in the charge sector of the many-body theory is io c q = t>F<?y 1 + 2V(aq)/irvF, 
whereas the spin degrees of freedom propagate with the free dispersion relation uo s q = vpq. 
We see that in the low-energy limit these degrees of freedom are completely separated as in 
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TL model. However, in the present case the charge excitation energy is not a linear function 
of q because the interaction behaves as V(aq) ~ j- In \ aq\ for \aq\ — > 0. Nevertheless, u) c q — > 
as 9 — ► 0. The spin degrees of freedom are not affected by the interaction since Viaq) only 
couples the charge density fluctuations. 

The model with the logarithmic divergence of V(aq) leads to non- analytic properties of 
thermal quantities. For example, the specific heat of this system at low temperatures is 
C = 7T + /3TlnT contrasting with the standard result in TL model where Ctl = JtlT^ 
Additionally, a single-particle density of states vanishes at the Fermi level as N(u) ~ u x ln y u 
with x and y being non-universal constants. Again in TL model Ntl{^) ~ cu^.i Also, as 
shown by Schultz,! the long range correlation functions have logarytmic corrections, and 
this might drive the system into a Wigner crystal. 

In conclusion, these results strongly suggest that the ground state and the low-energy 
properties of ID electrons in a quantum wire with the realistic form of the interaction are 
very different from those expected in the framework of the standard TL-type models. In 
particular, the actual form of the potential V(aq), and particularly, the screening effects due 
to the boundaries should modify transport properties in this system. It would, therefore, be 
very interesting to evaluate directly the conductivity. 
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FIGURES 

FIG. 1. Fourier transform of the electrostatic interaction, V{aq) as a function of aq for different 
ratios e\/e2 of the dielectric constants inside and outside the wire of the radius a. The inset shows 
V(aq) as a function of €2/^1 for \aq\ =0.1 in our model (||) (solid line) and in the model without 
the image charges (dashed line). 

FIG. 2. Exchange energy as a function of e±/e2 for different: a) radius a = 1,2,5,10 with 
r s = 1.0, b) RPA parameters r s = 0.5, 1, 2,4 (from bottom to top) with a = 1.0. Energy and length 
units are given in the units of the effective Rydberg and Bohr radius calculated with the effective 
mass m* inside the wire and the dielectric constant €2 outside the wire. The inset compares our 
model (Eq. |j) (solid lines) with the model without the image charges (dashed lines). 

FIG. 3. Correlation energy as a function of e\/e2 for different: a) radius a = 1,2,5,10 with 
r s = 1.0, b) RPA parameters r s = 0.5, 1, 2, 4 with a = 1.0. 

FIG. 4. Critical value of the RPA parameter r s as a function of ei/e2 for a = 1, 10. Above r c s 
the ferromagnetic phase is a stable ground state. Inset shows the behavior of r c s for small e\ < 1. 
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Fig.1, Byczuk & Dietl, PRB 
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Fig. 2a, Byczuk & Dietl, PRB 
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Fig. 3a, Byczuk & Dietl, PRB 
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Fig. 3b, Byczuk & Dietl, PRB 
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Fig. 4, Byczuk & Dietl, PRB 
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